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Spectral Analysis of the Stromlo-APM Survey II. 

Galaxy luminosity function and clustering by spectral type 
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ABSTRACT 

We study the luminosity function and clustering properties of subsamplcs of local 
galaxies selected from the Stromlo-APM survey by the rest-frame equivalent widths of 
their Ha and [O n] emission lines. The bj luminosity function of star-forming galaxies 
has a significantly steeper faint-end slope than that for quiescent galaxies: the majority 
of sub-L* galaxies are currently undergoing significant star formation. Emission line 
galaxies are less strongly clustered, both amongst themselves, and with the general 
galaxy population, than quiescent galaxies. Thus as well as being less luminous, star- 
forming galaxies also inhabit lower-density regions of the Universe than quiescent 
galaxies. 
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1 INTRODUCTION 

Important clues to the physics of galaxy formation and evo- 
lution may be obtained by studying the global properties, 
such as the luminosity function and correlation function, 
of quiescent versus star-forming galaxies. The most reliable 
tracer of the formation rate of massive, hot stars is the flux 
of the Ha emission line, directly related to the stellar UV 
(< 912 A) photoionizing flux (Kennicutt 1983). This line is 
frequently redshifted out of the observed spectral window, 
and so most deep galaxy surveys have instead used the [O n] 
3727A line as a measure of star-formation rate (Kennicutt 
1992). 

The luminosity function of galaxies subdivided by the 
presence or absence of the [O n] emission line has been cal- 
culated in the local Universe for the Las Campanas Red- 
shift Survey (LCRS, Lin et al. 1996a) and for the ESO Slice 
Project (ESP, Zucca et al. 1997). In both surveys it was 
found that the faint-end of the galaxy luminosity function is 
dominated by [O n] emitters, in other words that presently 
star-forming galaxies tend to be less luminous than quies- 
cent galaxies in both the bj (ESP) and Gunn-r (LCRS) 
bands. These results from [O ll]-selected samples are consis- 
tent with the recent luminosity function estimates from local 
samples of galaxies selected by morphological (eg. Marzke 
et al. 1998) and spectral (eg. Bromley et al. 1998, Folkes et 
al. 1999) type: early-type (elliptical and lenticular) galax- 



ies tend to be luminous, and late-type (spiral and irregular) 
galaxies faint. 

It is by now also well known (eg. Davis & Geller 1976, 
Giovanelli et al. 1986, Iovino et al. 1993, Loveday et al. 1995) 
that galaxies of early morphological type cluster together 
on small scales more strongly than late-type galaxies. Since 
emission- line galaxies (ELGs) tend to be of late Hubble type, 
we would expect ELGs to be more weakly clustered than 
non-ELGs, and indeed this has been observed by numer- 
ous authors (eg. Iovino et al. 1988, Salzer 1989, Rosenberg, 
Salzer & Moody 1994 and Lin et al. 1996b). 

In this paper we study the luminosity function and clus- 
tering for subsamples of the Stromlo-APM survey (Loveday 
et al. 1996) selected by Ha and [O n] emission-line equivalent 
widths. The Stromlo-APM survey is ideal for quantifying the 
statistical properties of emission-line versus quiescent galax- 
ies in the local universe since it contains a representative 
sample of different galaxy types and covers a large volume 
V « 1.38 x 10 6 fo~ 3 Mpc 3 . Since the red wavelength cover- 
age of Stromlo-APM spectra extends from 6300-7600 A we 
are able to detect the Ha (6562. 82A) line, when present, 
to a redshift z £ 0.16, i.e. beyond the maximum distance 
reached by the survey. Thus for the first time we are able 
to classify a large, representative sample of galaxies by the 
primary tracer of massive star formation, viz. the equivalent 
width of the Ha emission line. Measurement of the spectral 
properties of Stromlo-APM galaxies is discussed by Tresse et 
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al. (1999), hereafter referred to as Paper 1. The subsamples 
selected by their emission-line properties are described in 
The luminosity functions of the different samples are com- 
pared in ^ and in §^ we present clustering measurements. 
We summarize our results in Throughout, we assume 
a Hubble constant of Ho = lOO/i km s _1 Mpc _1 with h = 1 
and a deceleration parameter qo = 0.5. The exact cosmology 
assumed has little effect at redshifts z ^ 0.15. 



2 GALAXY SAMPLES 

Our sample of galaxies is taken from the Stromlo-APM red- 
shift survey which covers 4300 sq-deg of the south galactic 
cap and consists of 1797 galaxies brighter than bj = 17.15 
mag. The galaxies all have redshifts z < 0.145, and the mean 
is (z) — 0.051. A detailed description of the spectroscopic 
observations and the redshift catalog is published by Love- 
day et al. (1996). 

Of the 1797 galaxies originally published in the redshift 
survey, 82 have bj < 15. These bright galaxies are excluded 
from our analysis since they tend to be saturated on the 
Schmidt plates and hence have unreliable magnitudes. Of 
the remaining 1715 galaxies, 26 have a redshift taken from 
the literature, and for 7 we could not retrieve the spectra be- 
cause they were not observed with the Dual-Beam Spectro- 
graph (DBS) of the ANU 2.3-m telescope at Siding Spring. 
Also excluded were 6 blueshifted spectra, 3 with cz < 1000 
km s - , and 2 with too low signal-to-noise. 

The remaining 1671 spectra were flux-calibrated and 
had their spectral properties measured as described in Pa- 
per 1. Flux calibration of our spectra is accurate to ~ 10- 
20%, and so in the present paper we have restricted our 
analysis to galaxy samples selected by the equivalent widths 
(EWs) of their Ha and [O n] emission lines, which are in- 
sensitive to flux calibration errors. Note that since the res- 
olution of our spectra has FWHM = 5A, the Ha line can 
always be deblended from the [Nil] doublet. 

Of the 1671 measured galaxies, 11 were not part of our 
core statistical sample, either because they had an uncer- 
tain redshift or happened to lie in a part of the sky masked 
by "holes" around bright stars, etc. Of the remaining 1660 
galaxies, 82 could not have EW (Ha) measured as their red- 
shift places the Ha line in a small gap in the red part of 
the spectrum from 7000-7020A (Loveday et al. 1996). For 
an additional 57 spectra, Ha was seen in emission but could 
not be measured due to contamination by a sky line, or 
some other problem with the spectrum; [O n] lines could 
not be measured for similar reasons for 5 spectra. Note that 
lack of EW measurement, while correlated with redshift, is 
uncorrelated with galaxy morphology, and so we can reli- 
ably correct for missing EW measurements. We are thus left 
with a sample of 1521 galaxies which could be analysed by 
EW (Ha), and 1655 which could be analysed by EW ([O 
Histograms of log EW (Ha) and log EW([0n]) are plotted 
in Figures jj] and ^| respectively. 

We select galaxy subsamples using measured equivalent 
widths of the Ha and [O n] emission lines. The Ha line is 
the best tracer of massive star formation (Kennicutt 1983) 
but we also select samples using the equivalent width of the 
[O n] line, as this line allows us to compare with other sur- 
veys in which Ha is not always within the wavelength range 
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Figure 1. Histogram of log EW (Ha) for all galaxies and for mor- 
phologically selected subsamples as labeled. Note that galaxies 
with no detected Ha are not shown here. The median EW (Ha), 
including non-detections, is given for each sample after the mor- 
phological type. 



measured. The Ha line is detected with EW > 2A in 61% of 
galaxies. Of these emission-line galaxies, half have EW (Ha) 
> 15 A. Thus we form three subsamples of comparable size 
by dividing the sample at EW (Ha) of 2A and 15A. In the 
case of the [O n] line, 60% of galaxies have EW > 2A, and 
of these half have EW([On]) > 9.6A. The galaxy samples 
selected by Ha and [O n] equivalent widths are defined in 
Table [lj. 

Most galaxies in the Stromlo-APM survey have had a 
morphological type (elliptical, lenticular, spiral or irregular) 
assigned by visual inspection of the galaxy image (Loveday 
1996, Loveday et al. 1996). In Table III we give the numbers 
of galaxies of each morphological type in each spectroscopi- 
cally selected subsample. In Figures l] and ^ we also plot the 
distribution of equivalent widths for these morphologically- 
selected subsamples. The sample labeled "Unk" consists of 
galaxies to which no morphological classification was as- 
signed. We see that early-type galaxies dominate when Ha 
or [O n] emission is not detected and are underrepresented 
when emission lines are detected. Conversely, the number 
of irregular galaxies increases significantly in the spectro- 
scopic samples which show strongest star formation. Strong 
star formation is known to disrupt the regularity in the 
shape of a galaxy. In the deeper universe, the apparent in- 
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Figure 2. Histogram of log EW([Oll]) for all galaxies and for 
morphologically selected subsamples as labeled. Note that galax- 
ies with no detected [O n] are not shown here. The median 
EW ([O II]), including non-detections, is given for each sample af- 
ter the morphological type. 



Table 1. Spectroscopic subsamples and correlation with morpho- 
logical type. 



Sample EW (Ha) E SO Sp Irr Unk Total 

(a) H-low < 2 A 125 108 207 10 149 599 

(b) H-mid 2-15 A 8 16 340 18 81 463 

(c) H-high > 15 A 11 9 303 41 95 459 

Sample EW([On]) E SO Sp Irr Unk Total 

(d) O-low < 2 A 120 112 239 8 177 656 
(c) O-mid 2-9.6 A 19 24 344 19 97 503 
(f) O-high > 9.6 A 12 12 339 47 86 496 



galaxies have no morphological type assigned (the column 
marked "Unk" in Table [j]). The low median EW (Ha) and 
EW ([O n]) for these unclassified galaxies compared with the 
total sample suggests that many are in fact of early morpho- 
logical type. The spectral classification described in this sec- 
tion allows these galaxies to be assigned to their appropriate 
class in a quantitative way. 



3 THE GALAXY LUMINOSITY FUNCTION 

We estimate the b.j luminosity function (LF) for each galaxy 
subsample using maximum-likelihood, density-independent 
methods, so that our results are unbiased by galaxy cluster- 
ing. We use the Sandage, Tammann & Yahil (1979) paramet- 
ric maximum-likelihood estimator to fit a Schechter (1976) 
function, 

0(L)dL = 0*(-^) Q exp(^)dL. (1) 

We correct for random errors in our magnitudes by convolv- 
ing this luminosity function with a Gaussian with zero mean 
and rms a m = 0.30 (see Loveday et al. 1992, hereafter L92, 
for details). We also perform a non-parametric fit to each 
luminosity function using the stepwise maximum-likelihood 
estimator of Efstathiou, Ellis & Peterson (1988). This esti- 
mator calculates 4>{L) in a series of evenly-spaced magnitude 
bins and provides a reliable error estimate for each bin by 
inverting the information matrix, if-corrections are applied 
to each galaxy according to its morphological classification 
as E/S0: 4.14z, Sp: 2.252, Irr: 1.59z, Unk: 2.902. 

Before calculating the LF for each spectroscopic sub- 
sample defined in Table [l| we first checked that the galaxies 
omitted from this analysis, ie. those galaxies whose Ha or 
[O n] emission lines could not be measured, did not bias the 
LF measurement relative to the full Stromlo-APM survey. 
The LF estimates using all galaxies except the 194 with 
no Ha measurement available and all galaxies except the 
60 with no [O n] measurement were indeed both consistent 
with the full sample. 

Our estimates of the luminosity function for the 
EW (Ha) selected samples are shown in Figure |j| The inset 
to this Figure shows the likelihood contours for the best- 
fit Schechter parameters a and M* . The Schechter param- 
eters and their la errors (from the bounding box of the 
la error contours) are also listed in Table ^. Note that the 
estimates of a and M* are strongly correlated and so the 
errors quoted for a and M* in the Table are conservatively 
large. We see a trend of faintening M* and steepening a as 
EW (Ha) increases. There is a significantly greater contrast 
between the H-high and H-mid samples than between the 
H-mid and H-low samples, despite the rather similar dis- 
tribution of morphological types in the H-high and H-mid 
samples as compared with the H-low sample. This suggests 
that either there is not a simple one-to-one correlation be- 
tween optical morphology and EW(Ha), or that the larger 
fraction of Irr galaxies in the H-high sample are contributing 
to the steep faint-end slope for this sample. 

Luminosity function estimates of the EW([On]) se- 
lected samples and errors in the best-fit Schechter param- 
eters are shown in Figure ^. The la error contours for the 
O-low and O-mid samples overlap and the O-high sample 



crease in number of irregulars is also related to strong star 
formation (Brinchmann et al. 1998). Thus as expected we 
find a good correlation between morphological types and 
emission line equivalent widths. Since they can be measured 
objectively, spectroscopic properties of galaxies are a more 
reliable discriminator than visually assigned morphologi- 
cal types. Moreover, a significant fraction of Stromlo-APM 
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Table 2. Luminosity function parameters. 



Sample 


(V/Vmax) 


a 


M* 


n 


<t>* 


PL 


(a) H-low 


0.48 ±0.01 


-0.75 ± 0.28 


-19.63 ± 0.24 


10.1 ± 2.5 


4.5 ± 1.1 


5.9 ± 1.4 


(b) H-mid 


0.49 ±0.01 


-0.72 ± 0.29 


-19.28 ± 0.23 


11.0 ± 2.8 


5.4 ± 1.4 


5.1 ± 1.4 


(c) H-high 


0.54 ±0.01 


-1.28 ± 0.30 


-19.04 ±0.26 


48.0 ±15.9 


8.5 ±2.8 


8.8 ±2.9 


(d) O-low 


0.51 ±0.01 


-0.80 ±0.29 


-19.51 ± 0.22 


13.8 ± 3.2 


5.8 ± 1.3 


6.8 ± 1.7 


(e) O-mid 


0.49 ±0.01 


-0.36 ± 0.34 


-19.16 ± 0.22 


7.4 ±1.8 


5.8 ± 1.3 


4.9 ± 1.1 


(f) O-high 


0.51 ±0.01 


-1.49 ± 0.26 


-19.49 ± 0.30 


46.0 ±15.9 


3.7 ± 1.2 


7.9 ±2.8 



a is the faint-end slope and M* the characteristic bj magnitude of the best-fit Schechter 

function, n is the space density of galaxies in the range —22 < M < — 15 and <j>* is the 
normalisation of the Schechter luminosity function, both in units of 10~ 3 /i 3 Mpc~ 3 . is 
the luminosity density integrated over the same magnitude range, in units of 

lO 7 L h 3 Mpc- 3 . 




Figure 3. Estimates of the luminosity function for galaxies with no significant detected Ha emission (H-low: circles, solid line), with 
moderate Ha emission (H-mid: squares, dashed line) and with strong Ha emission (H-high: triangles, dot-dashed line). The symbols with 
error bars show the stepwise fit, the curves show the Schechter function fits. For clarity, data points representing fewer than five galaxies 
have been omitted from the plot. The inset shows 1 & la likelihood contours for the best-fit Schechter parameters. 



does not show a fainter M* than non-emission line galaxies. 
However, the LF for the O-high sample does have a signif- 
icantly steeper faint-end slope than that for galaxies with 
only weak or moderate [O n] emission. 

The fact that we see a systematic dimming of M* with 
emission-line EW for the Ha-selected sample but not for 
the [O n]-selected sample is probably due to the fact that 
EW (Ha) is a measure of the fraction of ionizing photons 
from OB stars over the flux from the old stellar population 
emitted in the rest-frame R band which forms the contin- 
uum at Ha, while EW ([O n]) is normalised by the flux from 
relatively young stars (mainly type A) . Thus EW (Ha) is 



more sensitive to the current star formation rate and hence 
blue luminosity enhancement than EW ([O n]). 

Note that the LF estimate for late-type galaxies pre- 
sented by L92 does not have such a steep faint-end slope as 
we find here for strong emission-line galaxies. In L92 we com- 
bined galaxies classified as spiral or irregular as "late type" , 
and so not all of them have strong emission lines. The faint- 
end slope for early-type galaxies (L92) was much shallower 
than that measured here for galaxies with no emission lines. 
At least part of this difference is due to a bias in the mor- 
phological type dependent LFs of L92 due to the tendency 
of unclassified galaxies in the Stromlo-APM survey to be of 
low luminosity (Marzke et al. 1994, Zucca et al. 1994). We 
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Figure 4. Estimates of the luminosity function for galaxies with no significant detected [On] emission (O-low: circles, solid line), with 
moderate [On] emission (O-mid: squares, dashed line) and with strong [On] emission (O-high: triangles, dot-dashed line). The symbols 
with error bars show the stepwise fit, the curves show the Schechter function fits. For clarity, data points representing fewer than five 
galaxies have been omitted from the plot. The inset shows 1 & 2<r likelihood contours for the best-fit Schechter parameters. 



avoid this bias with the spectroscopically selected samples 
analysed here. 

The normalisation <f>* of the fitted Schechter functions 
was estimated using a minimum variance estimate of the 
space density n of galaxies in each sample (Davis & Huchra 
1982, L92). We corrected our estimates of n, <j>* and luminos- 
ity density pL to allow for those galaxies excluded from each 
subsample. First, all subsamples were scaled by the factor 
1715/1660 to account for the 55 galaxies with no EW infor- 
mation available. Second, all Ha selected subsamples were 
scaled by 1660/1578 to account for the 82 galaxies whose 
Ha line, if present, would have fallen in the "red gap" (|0). 
Samples H-mid & H-high were scaled by an additional fac- 
tor 1578/1521 to allow for the 57 galaxies in which Ha was 
seen, but was not able to be measured. Finally, samples O- 
mid & O-high were scaled by 1660/1655 to allow for the 
five galaxies in which [O n] was seen but not measured. Our 
final estimates of n, <j>* and pL are given in Table |^. The 
uncertainty in mean density due to "cosmic variance" (L92 
equation 7) is « 6% for each sample. However, the errors 
in these quantities are dominated by the uncertainty in the 
shape of the LF, particularly by the value of the estimated 
characteristic magnitude M* . 

Using both Ha and [O n] equivalent widths as indicators 
of star formation activity, we find that galaxies currently un- 
dergoing significant bursts of star formation dominate the 
faint-end of the luminosity function, whereas more quies- 
cent galaxies dominate at the bright end. This is in agree- 
ment with the results of Lin et al. (1996a) and Zucca et al. 



(1997), but in disagreement with Salzer (1989), who finds no 
significant difference in the LF shapes of star-forming and 
quiescent galaxies. As pointed out by Schade & Ferguson 
(1994), Salzer's sample is biased against weak-lined ELGs 
at low- luminosity, and their reanalysis of his data correcting 
for this selection effect does find a steep faint-end slope for 
the LF of star-forming galaxies. 

The characteristic magnitude M* for the O-high sample 
is about 0.5 mag brighter than that for the H-high sample. 
This is probably due to a combination of several factors: 1) A 
large [O n] EW can come from a small [O n] flux and a very 
red continuum (ie. a small star formation rate and an old 
stellar population). 2) The correlation between estimated 
values of faint-end slope a and characteristic magnitude M* 
means that the steeper a of the O-high sample will push the 
estimated M* to brighter magnitudes. 3) The errors on M* 
are large (±0.3 mag), and so the H-high and O-high M* 
estimates disagree only at the 1-2 a level. 



4 GALAXY CLUSTERING 

In this section we measure the clustering properties of 
the galaxy subsamples. We measure the auto-correlation 
function of each sample in redshift space, and the cross- 
correlation function of each galaxy sample with all galaxy 
types in real space. For both estimates, we first verified that 
the 194 galaxies missing EW (Ha) measurement and the 
60 galaxies missing EW([Oll]) did not bias the measured 
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Figure 5. The solid histograms show the observed number- 
distance N(x) distribution for each galaxy subsample defined in 
Table [l]. The dashed lines show the expected distribution from 
our luminosity function fits. 



clustering relative to the complete sample. Those galaxies 
excluded because Ha fell in the "red gap" lie at redshifts 
z m 0.06-0.07. Nevertheless, omitting these galaxies did not 
significantly affect the measured clustering in real or redshift 
space. 



4.1 Redshift- Space Correlation Function 

We correct for boundary conditions and the survey selection 
function by populating the survey volume with a catalogue 
of ~ 18, 000 random points whose radial density matches 
that expected for each subsample. The number-distance dis- 
tributions for the six galaxy subsamples analysed here are 
shown in Figure ^. These plots also show the expected distri- 
butions inferred from the luminosity functions calculated in 
the previous section. We see that given the tendency for non- 
ELGs to be luminous and for ELGs to be faint, the ELGs 
are slightly overdense at large distances (a; J5, 200/i _1 Mpc) 
whereas there is an underdensity of non-ELGs at similar 
distances. This observation is reflected by the increasing 
(V/Vmax) with EW (Ha) seen in Table ^, and is probably 
due to evolution in emission line strength with redshift (eg. 
Broadhurst et al. 1992), occurring at redshifts as low as 
2 0.15. It is unlikely to be due to the changing pro- 



Figure 6. Estimates of the redshift-space correlation function 
for the galaxy samples given in Table |l| Error bars show the 
rms variance from dividing the survey into 4 distinct zones. The 
dashed line shows the best- fit power-law over the range 1.5—30 
h~ 1 Mpc with the index held fixed at 7 S = 1.47. The dotted line 
shows £ (s) estimated from the full Stromlo-APM sample (Loveday 
ct al. 1995). 



jected size of the spectrograph slit at different redshifts as we 
demonstrated in Paper 1. We checked that these discrepan- 
cies between observed and expected N(x) distributions did 
not bias our estimates of £(s) by also generating a random 
distribution according to a fourth-order polynomial fit to 
the observed radial density of each subsample. Clustering 
estimates using this random distribution gave results con- 
sistent with a random distribution generated according to 
the predicted radial density. 

The auto-correlation function of each sample in redshift 
space is measured using the estimator 



!+£(*) = 



■Wgg(s)w rr (s) 

[w gr (s)] 2 



(2) 



Hamilton (1993). Here w gg (s), w gr (s) and w rr (s) are the 
summed products of the weights of galaxy-galaxy, galaxy- 
random and random-random pairs respectively at separation 
s. We use the minimum- variance pair weighting given by 
equation 1 of Loveday et al. (1995), and the reader is referred 
to that paper for further details. Errors are estimated by 
dividing the survey into four zones of roughly equal area 
and calculating the variance in £(s) from zone-to-zone. 
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Table 3. Correlation function parameters. 



Sample 


so 


7r 


^0 


(a) H-low 


8.7 ±0.5 


1.78 ±0.08 


6.0 ± 1.4 


(b) H-mid 


5.5 ± 0.7 


1.60 ±0.13 


5.2 ± 2.0 


(c) H-high 


4.6 ± 0.9 


1.87 ±0.16 


2.9 ± 1.9 


(d) O-low 


8.6 ± 1.1 


1.79 ±0.07 


6.2 ± 1.8 


(e) O-mid 


4.9 ± 0.6 


1.64 ±0.05 


4.7 ±0.8 


(f) O-high 


4.1 ± 0.9 


1.78 ±0.15 


2.9 ± 0.7 



so is the correlation length measured in redshift space over the 
range 1.5-30 /i _1 Mpc with the power-law index held fixed at 
7 S = 1.47. 7r and ro are the real-space power-law parameters 

over 0.2—20 h~ x Mpc determined from cros s-correlation with the 
2d APM survey (jO). 



Estimates of £(s) are shown in Figure |(| A power-law 
£(s) = (s/so)~ 7s was fitted over the range 1.5-30 ft, - Mpc. 
For each subsample the estimated power-law slope 7 S was 
formally consistent with 7 S = 1.47, measured for the whole 
Stromlo-APM sample (Loveday et al. 1995) . Since estimates 
of the index 7 S and correlation length so are strongly cor- 
related, we determined the best fit so to each subsample, 
keeping the power-law index fixed at 7 S = 1.47. The results 
of these fits are shown by the dashed lines in Figure |^ and 
the best-fit values of so with la uncertainties (determined 
from fitting to each zone separately) are shown in Table |^. 

We see that the correlation length so becomes signif- 
icantly smaller in more actively star-forming galaxies, as 
traced by both EW (Ha) and EWQOll]). This result is 
in agreement with the power-spectrum analysis of the Las 
Campanas Redshift Survey by Lin et al. (1996b) who find 
that the clustering amplitude of ELGs is only about 70% 
that of the full LCRS sample. These results are also con- 
sistent with those of Rosenberg et al. (1994), Iovino et al. 
(1988) and Salzer (1989), all of whom find that ELGs are less 
strongly clustered than quiescent galaxies. Galaxies with no 
detected Ha (H-low) or [O n] (O-low) emission have a cor- 
relation length about twice that of ELG galaxies (H-high 
and O-high samples). This is larger than the difference in 
clustering amplitude determined by Lin et al. (1996b) from 
the LCRS, presumably because we have subdivided galaxies 
into three EW bins compared to their two EW bins. 



4.2 Real-Space Correlation Function 

The estimate of £(s) described above is affected by redshift 
space distortions. On small scales, random, thermal motions 
tend to decrease galaxy clustering, whereas on large scales, 
galaxy streaming motions tend to enhance £(s). In order 
to avoid the effects of galaxy peculiar velocities, we have 
calculated the projected cross-correlation function H(cr) of 
each galaxy subsample with all galaxies in the APM survey 
to a magnitude limit of bj = 17.15. We then invert this 
projected correlation function to obtain the real space cross- 
correlation function £(r) of each subsample with the full 
galaxy sample. This method of estimating £(r) is described 
by Saunders et al. (1992) and by Loveday et al. (1995). 

The large number of galaxy pairs used by this estimator 
allows us to fit a power-law to the measured cross-correlation 
function over the range of separations 0.2-20 ft _1 Mpc and to 
fit both the power-law index y r and the correlation length ro. 
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Figure 7. Estimates of the real-space correlation function for the 
galaxy samples given in Table |l| Error bars show the rms vari- 
ance from dividing the survey into 4 distinct zones. The dashed 
line shows the best-fit power-law over the range 0.2—20 /i -1 Mpc. 
The dotted line shows £(r) estimated from the full Stromlo-APM 
sample (Loveday et al. 1995). 



Our estimates of £(r) are plotted in Figure (?] and our best-fit 
power-laws are tabulated in Table |^. As in redshift-space, 
we see that strong emission-line galaxies are more weakly 
clustered than their quiescent counterparts by a factor of 
about two. 

The real space clustering measured for non-ELGs is very 
close to that measured for early-type (E + SO) galaxies, and 
the clustering of late-type (Sp + Irr) galaxies lies between 
that of the moderate and high EW galaxies (cf. Loveday 
et al. 1995). Given the strong correlation between morpho- 
logical type and presence of emission lines (Table [j]) this 
result is not unexpected. The power-law slopes are consis- 
tent (jr = 1-8 ± 0.1) between the H-low, H-high, O-low and 
O-high samples. For the moderate EW galaxies (H-mid and 
O-mid samples) we find shallower slopes (7,- = 1.6 ± 0.1). 
This is only a marginally significant (1-2 a) effect, but may 
indicate a deficit of moderately star-forming galaxies princi- 
pally in the cores of high density regions, whereas strongly 
star forming galaxies appear to more generally avoid over- 
dense regions. 
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5 CONCLUSIONS 

We have presented the first analysis of the luminosity 
function and spatial clustering for representative and well- 
defined local samples of galaxies selected by EW(Ha), the 
most direct tracer of star-formation. We have also selected 
galaxies by EW([Oll]), and find broadly consistent results 
between the two tracers of star formation, which is expected 
from their close relation (Kennicutt 1992, Paper 1). The 
observed trend for M* to fainten systematically with in- 
creasing EW (Ha), contrasted with the roughly constant M* 
with varying EW([Oll]), is probably due to EW (Ha) be- 
ing a more reliable indicator of star formation rate than 
EWQOll]). 

Star-forming galaxies are likely to be significantly 
fainter than their quiescent counterparts. The faint-end 
(M <; M*) of the luminosity function is dominated by ELGs 
and thus the majority of local dwarf galaxies are currently 
undergoing star formation. 

Star-forming galaxies are more weakly clustered, both 
amongst themselves, and with the general galaxy popula- 
tion, than quiescent galaxies. This weaker clustering is ob- 
servable on scales from 0.1-10 /i -1 Mpc. We thus confirm 
that star-forming galaxies are preferentially found today in 
low-density environments. 

A possible explanation for these observations is that lu- 
minous galaxies in high-density regions have already formed 
all their stars by today, while less luminous galaxies in low- 
density regions are still undergoing star formation. It is not 
clear what might be triggering the star formation in these 
galaxies today. While interactions certainly enhance the rate 
of star formation in some disk galaxies, interactions with lu- 
minous companions can only account for a small fraction of 
the total star formation in disk galaxies today (Kennicutt 
et al. 1987). Telles & Maddox (1999) have investigated the 
environments of H II galaxies by cross-correlating a sample 
of Hn galaxies with APM galaxies as faint as bj = 20.5. 
They find no excess of companions with Hi mass J> 10 s Mq 
near H II galaxies, thus arguing that star formation in most 
Hn galaxies is unlikely to be induced by even a low- mass 
companion. 

Our results are entirely consistent with the hierarchical 
picture of galaxy formation. In this picture, today's lumi- 
nous spheroidal galaxies formed from past mergers of galac- 
tic sub-units in high density regions, and produced all of 
their stars in a merger induced burst, or series of bursts, 
over a relatively short timescale. The majority of present- 
day dwarf, star-forming galaxies in lower density regions 
may correspond to unmerged systems formed at lower peaks 
in the primordial density field (eg. Bardeen et al. 1986) and 
whose star formation is still taking place. Of course, the full 
picture of galaxy formation is likely to be significantly more 
complicated than this simple sketch, and numerous physical 
effects such as depletion of star-forming material and other 
feedback mechanisms are likely to play an important role. 
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